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Magnetic twist of the active region has been measured over a decade using 
photospheric vector field data, chromospheric H Q data, and coronal loop data. 
The twist and tilt of the active regions have been measured at the photospheric 
level with the vector magnetic field measurements. The active region NOAA 
10930 is a highly twisted emerging region. The same active region produced 
several flares and has been extensively observed by Hinode. In this paper, we 
will show the evolution of twist and tilt in this active region leading up to the 
two X-class flares. We find that the twist initially increases with time for a few 
days with a simultaneous decrease in the tilt until before the X3.4 class flare 
on December 13, 2006. The total twist acquired by the active region is larger 
than one complete winding before the X3.4 class flare and it decreases in later 
part of observations. The injected helicity into the corona is negative and it is 
m excess of 10 43 Mx 2 before the flares. 
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1. Introduction 

It is generally believed that the magnetic fields on the sun is generated 
underneath the convection zone and they rise due buoyancy into the pho- 
tosphere. Active region magnetic fields emerge as shaped fluxtubes at the 
photosphere. These emerged magnetic fields always exhibit some amount 
of twist. Twist is shearing of the field lines about its axis. Leka et al. 1 have 
found that the twisted emerging flux carries a current along with it. The 
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twist in the active region can originate due to several reasons. In the sub- 
photosphere, the held lines could twist due to the magnetic forces that are 
internal to fields. The generated magnetic fields can be twisted during the 
dynamo processes. 2 The generated field while rising through the convection 
zone is influenced by the convective buffeting that can introduce the twist 
in the field lines. 3 In the photosphere, the differential rotation of the sun 
can also introduce the twist, though in small amount. 4 The shearing motion 
of the magnetic footoints can also twist the field lines at the photosphere 
after the emergence. 

The parameter a appear in the force-free field equation VxJ] = aB is 
considered as a proxy for the twist. 5 Active region twist exhibits a hemi- 
spheric helicity rule, that is left/right handed twist in the northern and 
southern hemisphere. 6 ' 7 However, many researchers now believe that the 
hemispheric helicity rule has to be re-investigated due to inconsistencies 
reported for different phases of solar cycles. 8-10 

Tilt is the angle made by the line joining the centroids of two opposite 
polarity with the solar equator. Many researchers considered the tilt as a 
proxy for the writhe of a fluxtube. The writhe is the spatial deformation 
of the flux tube around its axis. The tilt of the active region also follows 
the hemispheric rule that is positive in the northern and negative in the 
southern hemisphere, following the Hale-Nicholson law. 11 The tilt also de- 
pends on the latitude and is known as Joy's law. 12 The twist and writhe is 
interchangeable and their sum is a conserved quantity. 13 

The twist parameter is an interesting quantity in emerging active re- 
gion. 14 Pevtsov et al. 14 have observed that the a increases as the active 
region emerges and reaches a constant value after a day of the flux emer- 
gence. As the flux emerges into the photosphere, the coronal helicity in- 
creases. The accumulated helicity in the corona plays an important role in 
the energetic events such as flares and coronal mass ejections. Nandy et 
al. 15,16 have studied the relaxation of the twist towards constant a force- 
free fields after flares. Nindos and Andrews 17 relate the coronal magnetic 
field line twist to eruptive and non-eruptive events. 

The active region NOAA 10930 is well studied for twist, 18 flares 19 etc. 
In this paper, we would like to show how the twist and tilt of the active 
region varies over the time. Also, we would like to investigate, if there is 
any changes in the twist during the flare or CME. We used 

Qmean &S & 

parameter to study the twist in the active region. In the next section, we 
describe the vector magnetogram data utilized in this study and in Section 
3, we present the results on how the twist and tilt evolves with time. In 
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Section 4, we present our conclusions. 



2. Observational Data and Its Analysis 

The spectro-polarimctcr is one of the back-end instruments in solar opti- 
cal telescope (SOT 20 - 22 ) on board HINODE 23 satellite that makes spec- 
tropolarimetric measurements at a spatial sampling of 0.3" in fast mapping 
mode and 0.16" in normal mode. The Stokes I, Q, U and V spectra in Fe 
I 6301.5 A and 6302.5 A lines are obtained. It makes a map of the active 
region by spatial scanning of the slit positions. The spatial resolution along 
the slit direction is 0.295" and in the scanning direction is 0.317"/pixel. 
The obtained data set consists of total 36 vector magnetograms and is a 
mixer of both the modes. The obtained Stokes signals are calibrated using 
the standard solar software pipeline for the spectropolarimetry. The com- 
plete information on the vector magnetic fields is obtained by inverting the 
Stokes vector using the Milne Eddington inversion. 24-26 The ambiguity in 
the transverse field is resolved based on the minimum energy algorithm 
developed by Metcalf 27 and Leka. 28 The transformation of the magnetic 
field vectors to the heliographic coordinates has also been done. 29 The er- 
ror in vertical (B z ) and transverse field (B t ) measurement is 8 G and 23 G, 
respectively. 

Several methods have been developed to measure the twist in the ac- 
tive region magnetic fields. Almost all methods are based on the force-free 
parameter a. The a^est was introduced by Pevtsov. 6 This is the value of 
alpha for which the computed transverse field is best matched with the 
observed transverse field for whole active region and is called as ab est ■ The 
other methods to measure the overall twist in the active region magnetic 
fields are (a) average value of a computed for overall active region magnetic 
field, a mean 30 (b) global alpha, a g and (c) peak value of a in the active 
region, a pea k- 31 Apart from these, there is one more method that has been 
introduced in the literature known as signed shear angle (SSA). 10 

In this paper, we use the method of a mean to measure the twist in the 
active region NO A A 10930. We express the total twist in the emerging 
bipolc region as a product of a magnetic coronal loop length (1) and the 
winding rate (q). If one takes the coronal loop as a semicircle of length '1' 
then the circumference of the semicircle is ird/2. Where d is the distance 
between the flux weighted centroids of the two poles. The winding rate q 
is related to the twist parameter as a mean /2. Then the total twist is the 
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product of lq 31 ' 32 and is given by, 

_ ~2 2 ^ ' 

The tilt angle is computed by measuring the tilt of the line joining 
the flux weighted centroids of positive and negative poles in the counter- 
clockwise direction with respect to the equator. The east side of the line 
joining the equator is zero and the west side is 180° with values increasing 
in the counter-clockwise direction as a measure of tilt in the active region 
NOAA 10930. The pixels with vertical field strength values larger than 
100 G is utilized in computing the centroid. This is to avoid the small scale 
features and noisy pixels. 

The alpha maps have been obtained by first computing the vertical 
current density as 

'■(*.»> 

once the vertical current density is estimated then the vertical compo- 
nent of the twist parameter is computed through the relation between a z , 
J z , and B z as, 

a z (x,y) = (3) 
B z 

where, B z is the vertical component of the magnetic field in heliographic 
co-ordinate system and /iq is the permeability of free space and its value 
is 47rxl0~ 3 G mA" 1 . To avoid the pixels with zero B z , that always occur 
near the neutral line, we have used a threshold of 100 G for B z . The a z 
map is useful for finding the evolution of twist in space and time. 

It is found that the measured twist (a z ) at the photosphere closely 
agrees with the coronal field line twist (a corona )- 5 In order to compute the 
coronal hclicity we followed the results obtained by Berger 33 and imple- 
mented to the observed magnetograms as done by Demoulin et al. 34,35 The 
relative magnetic helicity, using a constant alpha force free method in the 
linearized form is given by, 



H r = a V V ' "7^' (4) 

where, B nxiTly is the Fourier amplitude of the longitudinal component 
of the magnetic field, k x — 2im x /L, k y — 2im y /L, and L is the length of 
the computational window in the horizontal axis. The original form of the 
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Fig. 1 . The vector magnetograms are shown for each day. The background black (white) 
color represents the S (N) polarity regions of the vertical magnetic field strength (B z ). 
The arrows represents the transverse field vectors overlaid upon B z . The size of the 
arrow at the bottom of each figure indicates the magnitude of the transverse field and 
the vertical bar on the top right figure shows the magnitude of B z expressed in G. 
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equation, its limitations and its simplification to the above form can be 
found in Demoulin et al. 34 In estimating the coronal helicity we used the 
ctmean as the value for a in the above equation. 

3. Results 

The active region NOAA 10930 appeared in the declining phase of the solar 
cycle 23 rd in southern hemisphere at a latitude of 5° . The active region had 
two polarities when it appeared on the east limb on December 5, 2006 and 
they were almost aligned perpendicular to the East- West direction. It still 
survived when it crossed the west limb on December 17. 



3.1. Magnetic Flux 

Figure 1 show the sample images of transverse field vectors overlaid upon B z 
for each day of observations. In the vector field maps (as shown in Figure 
1) it is clear that field lines are twisted. Both the polarities are located 
close to each other and boundary between the two is highly sheared. The 
N-polarity region was small on December 09, 2006 and it evolved over a 
period of 7 days. 
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Fig. 2. The temporal evolution of magnetic flux is shown as a function of time. Top: 
N-polarity, Bottom: S-polarity. 



Figure 2 shows the temporal evolution of flux over a period of 7 days. 
From the plot it is clear that there is a large imbalance in flux in the N and S- 
polarities. A vertical dashed lines in the plot represents the time of X3.4 and 



December 2, 2010 1:38 



WSPC - Proceedings Trim Size: 9in x 6in ms 



7 

XI. 5 class flares that occurred on Dec 13, 2006 at 02:30 UT and Dec 14, 2006 
at 21:07 UT respectively. The flux in the N-polarity region continuously 
increased until the X3.4 class flare occurred. It remained almost constant 
for sometimes and decreased subsequently. 

3.2. Twist in Active Region 

Twist is not distributed uniformly over the sunspot. In order to see where 
exactly the contribution to the overall twist is coming from, we have made 
the a z maps. In a z maps, it is easy to see that there is a mixed a in the 
penumbra as seen by Su ct al. 10,18 There is also salt-pepper like mixed 
alpha in the umbra. However, a dominant negative sign of a is present in 
the emerging region, the N-polarity sunspot and also in a small portion of 
S-polarity sunspot. These dominant negative sign of a in S and N-polarity 
region is shown by arrows in top- left of Figure 3. This dominant negative 
sign of twist can also be seen in all the subsequent maps made at different 
times in both the N and S-polarities. In the N-polarity region, the area of 
the dominant negative sign in a increases with time till Dec 12 and later it 
decreases in size. 

Figure 4(left) shows the plot of total twist (*) and tilt (x) expressed in 
terms of winding number. The total twist is estimated using the equation 
(1), and the number of windings is estimated as T/2-k. From the plot it is 
clear that the total twist in the active region increases with time till the 
mid of Dec 12, 2006 and after that it starts to decrease with time. The 
active region has got more than one complete winding on Dec 11 which 
stays there till the mid of Dec 12, 2006. The dashed vertical lines in the 
plot indicate the time of X3.4 and XI. 5 class flares. Before the X3.4 class 
flare the active region acquired about 0.8 turns and it reduced to 0.5 turns 
of the complete winding after the flare. Similarly, before the XI. 5 class 
flare the active region has got close to one complete turn, and it decreased 
to 0.4 after the flare. We have estimated the error in measuring the total 
twist by varying the thresholds in Bz from 50 to 250 G in steps of 50 G. We 
estimated the mean value of total twist for those pixels whose value exceeds 
the thresholds in B z . With this method we estimated the error in measuring 
the total twist and it is found to be less than 5%. This small variations in 
total alpha will not affect our results. Apart from these, we also examine 
the temporal variations in total twist using the another method called a g 
and found a similar variations in the total twist as a function of time. 

The x symbol in Figure 4(left) indicates the tilt of the active region with 
respect to the E-W direction. It started to decrease from the beginning 
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of the observation till the beginning of the X3.4 class flare. It attained 
almost constant value after the X3.4 class flare until the end of December 
15, 2006. The tilt and twist behaved exactly opposite from December 9, 
2006 till December 12, 2006. However, they differ in their magnitudes. This 
behavior is similar to the results obtained by Holder et al 36 and Nandy. 37 
However, in the present case as the total twist increased with time, the tilt 
decreased. 

It has been observed that the emerging N-polarity was rotating in the 
anti-clockwise direction. 38 The other S-polarity did not display any rota- 
tion. In order to see how much the twist appeared in the N-polarity we 
have independently estimated the twist for the N-region. Figure 4(right) 
shows a plot of total winding of N-polarity over a period of 7 days. Simi- 
lar to the total twist of whole active region, the N-polarity also exhibited 
the increase in total twist starting from December 9, 2006 till mid day of 
December 12, 2006. It accumulated more than one complete winding just 
before the X3.4 class flare. The winding reduced after the X3.4 class flare, 
however it was still larger than one complete winding. Later, the total twist 
decreased and before the XI. 5 class flare it was just 3/4th of one complete 
winding. Again, after XI. 5 class flare there is small amount of decrease in 
total twist. It should be noted here that this is a mean value of the a z 
taken from distribution of a z . This means that there is a value of a z larger 
than the mean value and lower than the mean value. It also indicates that a 
significant number of pixels show their total twist larger than one complete 
winding. 
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Fig. 5. A Plot of evolution of relative magnetic helicity as a function of time is shown 
here. 

3.3. Relative magnetic helicity 

Due to the emergence and rotation of sunspots, there is a large amount of 
helicity accumulated in the solar corona. There are several ways to estimate 
the relative magnetic helicity in the solar corona. Since the cadence of vector 
magnetic field data is small we cannot compute the velocity vectors required 
to estimate the relative helicity injection. Hence, we have used the eq (4) 
to estimate the coronal magnetic helicity. We have used the a mean as a in 
the eq(4) . Figure 5 shows a relative magnetic helicity plotted over a period 
of 7 days. The magnetic helicity curve follows the total twist curve. The 
injected helicity is very large. Before the X3.4 class flare it was 7xl0 43 Mx 2 
which decreased to 4xl0 43 Mx 2 after the X3.4 class flare. 

4. Conclusions 

An increase in the flux and the total twist of the active region during 
first half of the observation suggests that a new flux region is emerging. 
This confirms that the emerging region carries current and twist generated 
below the photosphere along with it as observed by Leka et al. 39 The twist 
increased in the active region until the on-set of X3.4 class flare and showed 
a decreasing trend thereafter. A similar trend of increase and decrease in 
twist was noticed before and after the XI. 5 class flare too. This suggests 
that the twist in the active region slowly relaxing towards the lower value. 15 
Leamon et al. 32 suggest that the coronal loops become unstable and 
erupt only when their total twist exceeds the critical twist, 2ir. However, 
in their study they do not find any active region exceeding this critical 
twist. Many active regions exhibit their total twist of about 2n/3 and some 
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of them even did not exceed 2tt/6 value of one complete winding. In our 
observations, we find that the total twist exceeded one complete winding on 
December 11 and 12, 2006. However, it decreased just before the X3.4 class 
flare reaching 0.8 times one complete winding and before the XI. 3 class 
flare it is close to one complete winding. This result was obtained when we 
take the active region as a whole. On the other hand, when we computed 
the total twist in the N-polarity region alone, it exhibited larger than one 
complete winding and it was about 1.3 times the 1-complete winding before 
the X3.4 class flare. It decreased to 1.1 turn after the X3.4 flare. However, 
before the XI. 3 class flare the total twist was less than 37r/4 of 1-complete 
turn. This is only for the mean value of the total twist in the active region. 
However, the active region has a distribution of local twist (a) and many 
regions show larger twist than the mean of this distribution. This means 
that this active region shows more than one complete winding in many 
locations of the sunspots, exceeding the critical limit of 2ir. We do not know 
at this stage whether this active region erupted due to kink instability or 
due to some other reason. We need a coronal observations to find such kink 
instability. At present, from this study it appears that there is a decrease in 
the total twist after the flare. In order to firmly establish the role of critical 
twist in the eruptive X and M class flares, we need to observe many such 
events. 
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